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ABSTRACT 

In  work  are  received  the  oretical  results  on  blow-up  for  a  class  of  problems  with  additional  linear  and 
nonlinear  diffusion  and  it  would  be  interesting  to  understand  reaction-diffusion  profiles  for  these  problems  as  well. 
Method  of  characteristics  does  not  directly  apply  to  such  problems  due  to  the  diffusive  nature  of  the  dynamics.  Special 
research  techniques  of  the  nonlinear  parabolic  equations  which  allow  conducting  rather  detailed  research  blow  up  of 
solutions  of  a  heat  conduction  equation  with  a  source  were  developed. 

KEYWORDS:  Nonlinearity,  Diffusion  &  Reaction 


TRANS 

STELLAR 

•Journal  Publications  •  Research  Consultancy 


Received:  Apr  03,  2019;  Accepted:  Apr  23,  2019;  Published:  May  23,  2019;  Paper  Id.:  IJMPERDJUN20 19117 

1.  INTRODUCTION 

In  recent  years  began  to  pay  special  attention  to  the  unlimited  decisions  which  are  the  reason  of  existence 
of  an  energy  release,  chemical  reaction,  etc.  Such  decisions  arise  in  many  physical  processes  (for  example, 
combustion).  In  this  regard  in  recent  years  the  theory  of  blow  up  of  decisions  strongly  develops,  many  works  of 
A.A.  Samarsky,  S.P.  Kurdyumov,  A.P.  Mikhaylov,  V.A.  Galaktionov,  S.N.  Dimova  and  many  foreign  scientists 
are  devoted  to  this  question.  Blow  up  of  the  decision  were  called  by  decisions  with  the  aggravation  mode.  Special 
research  techniques  of  the  nonlinear  parabolic  equations  which  allow  to  conduct  rather  detailed  research  blow  up 
of  solutions  of  a  heat  conduction  equation  with  a  source  [1,2,  3]  were  developed. 

The  group  of  the  Bulgarian  scientists,  in  particular  S.N.  Dimova,  M.S.  Kaschiyev,  M.G.  Koleva  and  D.P. 
Vasileva  enter  received  important  scientific  result  recently.  In  this  model  at  a  HS  mode  approximation  with 
aggravation  to  the  S-mode  the  possibility  of  existence  of  waves  with  the  complex  structure  of  the  organization 
which  are  also  the  structures  attractors  described  by  invariant  and  group  decisions  is  open.  Half-width  of  these 
structures  waves  grows  over  time.  It  was  supposed  earlier  that  the  difficult  world  of  structures  corresponds  only  to 
the  LS  mode  with  the  reduced  half-width  at  the  prevailing  role  of  action  of  nonlinear  sources  in  comparison  with 
diffusion  processes.  And  in  S.N.  Dimova's  work  with  colleagues  mentioned  above  still  the  difficult  world  of  the 
solitonic  structures  waves  keeping  the  form  at  the  growing  half-width  [3]  is  open. 

In  article  [4]  are  proved  an  existence  theorem  and  uniqueness  of  the  generalized  solution  of  an  initial 
value  problem  for  the  equation 

ut = WOL  +v(u). 
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Boundary  value  problems  for  the  similar  equations  are  investigated,  for  example,  in  the  book  [5]  where  it  is 
possible  to  find  further  references. 

Presence  at  an  evolutionary  task  of  the  unlimited  decision  means  its  global  insolubility  on  time  (a  non-existence  of 
a  global  solution).  In  this  case  at  the  expense  of  an  intensive  energy  release  process  of  combustion  can  happen  in  the 
so-called  mode  to  aggravation.  In  other  words,  the  initial  value  problem  has  no  global  solution  on  time  and  in  some  instant 

t  —  T0  <  +°°  aggravation  moment),  amplitude  of  the  solution  becomes  infinitely  big: 

sup  u(t, x )  — »  +°°,  t  — >  T0 . 

X £Rn 

At  the  same  time  up  to  the  aggravation  moment  at  all  0  <  t  <  T0  non-negative  continuous  solution  of  a  task  is 
limited  and  also  generalized.  It  cannot  have  all  derivants,  however  function  |Vw(f,x)|  is  continuous  everywhere  in  the 

area  (0,  T0  )x  R  N  .  From  the  physical  point  of  view  it  means  the  continuity  of  a  heat  flux  equal  W  =—  |x|m|Vw|"  V  m  . 

The  nonlinear  effect  of  final  rate  of  propagation  of  temperature  indignations  was  for  the  first  time  found  by 
Zeldovich  Ya.B.  and  Kompaneets  A.  S.  [6],  and  effect  of  spatial  localization  of  shift  indignations  -  Martinson  L.K.  and 
Pavlov  K.B.  [7].  On  examples  of  exact  self-similar  decisions  in  [8]  it  was  shown  that  shift  indignations  can  extend  with  a 
final  speed  unlike  the  incompressible  Newtonian  ( n-1 )  and  pseudoplastic  fluids  (/;</),  in  which  the  speed  of  distribution 
of  shift  indignations  is  infinite. 

Kalashnikov  A.S.  [9]  established  the  phenomenon  of  full  cooling  for  final  time  in  case  of  strong  absorption 
(0  <  /?  <  1)  . 

In  [10]  studied  the  following  Cauchy  problem: 

(!) 

mQT=RN  x(0,T),N>l, 

u(x,0)  =  u0(x),  reR*,  w0(x)  >  0  for  constant  x  e  RN  ^2) 

In  the  case  of  fast  diffusion  0  <  m  +  p  —  1  <  2  estabilished  bilateral  evaluations  Lrj  . 

Here 

\_ 

u  =  u(t,x),  x  =  (xv...,xN),  |x|  =  (x12  +  ...  +  j4)2,  Om  =  («v...,hJ  (3) 

m  +  p  -  3  >  0,  p  —  1  >  0 

p(s^ ,  S  >0  -  decreasing,  continuous,  positive  function;  /?(())  =1 
Besides, 
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suppu0e  BRo  ={\x\  </?„},  zv*"<0°  (4) 

Typical  example  of  the  function  p  is 

yO  =  (l+  |jc|)  ' ,  /  >0-  (5) 

It  is  said  that  a  solution  of  equation  (1)  has  a  finite  speed  of  propagation  property  (FSP)  if  the  perturbation  of  the 
conditions  suppw(.,  tQ  )<  oo  at  some  point  in  time  t{]  >  0  follows  that  this  property  remains  for  all  instants  t  >  t{)  Otherwise 

say  that  the  carrier  of  the  solution  (1)  collapses  for  a  short  time  (ShTC).  Main  objective  of  this  work  is  clarification  of 
conditions  on  /?(|.ij)  >  at  which  ShTC  properties  for  the  solution  of  a  task  (1),  (2)  take  place.  Let's  note  that  ifp  =  1,  then, 

by  virtue  of  conditions  (3)  and  (4),  solutions  of  problem  (1),  (2)  have  the  property  finite  velocity  of  propagation. 
However  if,  for  example,  in  (5)  1  "it  is  too  big",  then  takes  place  ShTC. 

In  [11]  the  quasilinear  degenerating  parabolic  equation  with  the  non-uniform  density  investigated.  It  established 
that  depending  on  behavior  of  density  on  infinity  for  the  solution  of  an  initial  value  problem  take  place  or  property  of  final 
rate  of  propagation  of  indignations,  or  destruction  of  the  carrier  for  final  time.  The  following  second  mixed  problem  is 
considered: 


t-div{um~'\Du\r  2Dm)  =  0  in  Qt=CIx(0,T), 

(6) 

m~1\Du\P~2^:  =  0  on  ddx(0,T) 
dn 

(7) 

:(x,0)  =  u0(x),  xe£l, 

(8) 

where  fie  K*1  N  >2  -  unlimited  domain,  mgr  O  =  lol  -oo-  40-  noncompact  rather  smooth  border  Q-  n~ 

’  N  |  Iat  ’ 

external  single  normal  to3fix(0,7\),  r>O.Are  supposed  tom  +  p_3  (o,  p -1)0,  m  + p -2>  max 

u0(x)>  0  xe  Cl  and  u0e  L,  ,r  (<T)It  is  known  [12]  that  atm+  p—  3<0(6)  fall  into  to  the  equations  describing 
process  with  fast  diffusion.. 

In  article  [13]  the  initial  value  problem  of  rather  parabolic  equations  with  double  nonlinearity  of  the  following 
look  is  considered:  u  =  div[ua  |Z)a|m_1  Du  j  +  u p  >  whereO  <  m  +  a  <  1.  It  establishes  the  existence  and  non-existence  on 

the  whole  of  the  time  of  solving  this  problem  for  initial  data,  slowly  approaching  zero. 

The  following  Cauchy  problem  is  considered: 

ut  =  div[ua \Du\n  '  Du  j  +  up  (9) 


(x,t)e  Qt=Rnx{0,T),  T  >0,  N>1 
u(0,x)  =  u0(x),  xeR*1 


www.tjprc.ors 


SCOPUS  Indexed  Journal 


editor@tjprc.  org 


1098 


Muhamediyeva.  D.  K 


Here  it  assumed  that 


m  +  a <  1, m  >  0, m  +  a>  max 


|o,l  — ^±l|,p>land  u0(x)  -  it  is 


non-negative 


measurable  function  from  a  class  ]_,  .  At  in  +  OC  —  1  <  0  equality  (9)  without  source  falls  into  to  the  equations  of  fast 


l, loci. 


diffusion,  and  at  m  —  1  it  arises  in  a  plasma  physics  [12].  The  characteristic  feature  of  this  class  of  the  equations  is  the  lack 
of  property  of  final  rate  of  propagation  of  indignations.  It  is  well  known  [14]  that  even  in  the  simplest  cases  m  =  1,  (X  =  0 

there  are  solutions  (9)  which  beyond  all  bounds  grow  for  of  course  time,  that  is  exists  T  >  0 ,  such  that 


^°°’  t^T< oo. 


It  was  shown  that  if  p  <  p  —\  +  2  /  N ,  that  any  non-negative  solution  of  a  task  (9),  (10)  "blows  up"  for  final 

*  * 

time.  If  p  >  p  for  initial  functions,  enough  small  in  a  sense,  there  is  a  decision  in  general  on  time.  If  p  =  p  joins  in  a 

mode  case  with  aggravation. 

2.  PROBLEM  DEFINITION 


In  the  present  work  questions  of  global  resolvability  of  an  initial  value  problem  with  double  nonlinearity  and 
qualitative  properties  of  the  solution  of  a  task  on  the  basis  of  the  self-similar  analysis  are  investigated.  These  qualitative 
properties  of  the  considered  task  it  carried  out  because  of  a  research  of  qualitative  properties  of  the  self-similar  equation. 


Let's  consider  in  the  domain  Q  =  {(t,x) :  0<t,xe  R}  quasilinear  equation  of  reaction  diffusion  with  double 
nonlinearity 


du 

Ht 


d_ 

dx 


Du" 


p-2  ^ 

an 

dx 


-  ku  (l  —  J,  u  |f=0=  u0(x )  >  0 ,  xe  R 


(11) 


which  describes  reaction  diffusion  process  which  diffusion  coefficients  are  equal  D  li" 


where 


in,  p,  P  -  positive  real  numbers,  li  =  u(t ,  x)  >  0  -  required  solution. 


Let's  note  that  the  main  feature  in  studying  of  nonlinear  properties  is  receiving  various  types  of  estimates  of 
solutions,  and  then  on  their  basis  numerical  model  operation  of  a  task.  In  it  the  larger  role  played  by  the  self-similar  and 
approximate  and  self-similar  approaches  that  are  widely  presented  in  particular  in  A. A.  Samarsky,  S.P.  Kurdyumov,  V.A. 
Galaktionov,  A.P.  Mikhaylov's  works  [16]. 

For  creation  of  the  self-similar  equation,  the  algorithm  of  nonlinear  splitting  is  used  [17]. 

Let's  note  that  replacement  in  (11)  (11)  u(t,  x)  —  ek’v(T(t),  x)  will  lead  (11)  to  a  look: 


3v 

3r 


d_ 

dx 


Dv" 


P~2dv 

dx 


y 


(12) 


v|f=0  =v0(x)  =  u0(x). 
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Taking  into  account  (3)  equation  (2)  will  correspond  in  the  following  look: 
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dv_d_ 
dr  dx 


Dv" 


3v 

dx 


-k(rbvp+\ 

(j3-(m+  p)-  3) 


(13) 


where 


k1=k((m  +  p-3)kf  b  m  +  p-  3 


For  the  purpose  of  receiving  the  self-similar  equation  we  will  apply  a  splitting  algorithm  to  the  equation  (13) 
according  to  which  we  will  find  at  first  the  solution  of  an  ordinary  differential  equation  — -  =  —k{UbV^+' , 

dr  1 


r  Bk 

Then  we  have  v (z)  =  c(z  +  TQ)  y ,  T0  >0 , where c  =  H  1 

b  + 1 


i 

7  b  + 1 
,  y  = 


P 


Then  the  solution  of  the  equation  (13)  looked  for  in  a  look  v{t,  X )  =  v(l)w(T,  x\  a  function  T  —  T(t )  is  chosen 


r(t)  =  jv(m+p-3)(t)dt  = 


,l-[y(m+p-3)] 


Inf, 


/ (1  - /(m  +  p- 3)),  if  1  - y{m  +  p  —  3)^0, 

if  1  -  ^(m  +  p  -  3)  =  0, 
if  m  +  p  =  3  . 


Then  for  w(r,  x)  let's  receive  the  equation 


dw  _  d 
dr  dx 

where 


Dw"' 


dw 


dx 


P~ 2  ->  ^ 

aw 
dx  j 


+  y/{vi-wp+l)  ■ 


(14) 


1 


y/-\  (1  -[y(m  +  p-3)])r 


if  1  -  [/(m + p  —  3)  >  0, 


rc-<r<’»+P-3))'  jf  l-[r(m  +  p-3)  =  0. 

That  shows  an  invariance  of  the  transformation  given  above.  Let  us  consider  the  self-similar  solution 

w(x,x)  =  /(£),  £,  =  |x|  /  x1/p  .  (16) 

for  equation  (14). 

Then  substituting  (16)  in  (14)  in  the  case  1  —  [/(m+  p  —  3)  >  0,  regarding  to  / (^)it  is  easy  to  receive  the 


(15) 


self-similar  equation: 


d 


L(f)  =  ^-(fm-' 


df 


d£ 


p-2 


^  +  -^T  +  M(f-f/l+l)  =  0,p=  1 - — 

dg  p  dg  \-[/{m  + p -3)] 


(17) 
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Let's  be  engaged  in  creation  of  the  upper  solution  for  the  equation  (11).  If /3  —  [3  —  (p  +  m)\/  (p  —  1)  ^ 

r  /  _  a  /  _|_  £Y\  n 

equation  (17)  has  the  exact  solution  of  a  look  ■'-  ’  '  ~  ^  '+ 

where  n  =  (p  —  l)  /  (p  +  m  —  3) ,  y=pl(p- 1),  {b)+  =  max(0,Z?)  [19]. 

Properties  of  solutions  of  equation  (17)  in  the  case  p=2,  ni=  /investigated  in  detail  in  [3]. 

Theorem  1 

Let  u(0,  X)  <  Z+  (0,  x),  XE  R.  Then  for  the  solution  of  problem  (11)  in  the  domain  Q  the  estimate 
u(t,x)  <  Z±(t,x)  =  (T  +t )  7  f±(  ^  =  |x|  /  XUp  Here  f±{£)  function  defined  above. 

3.  RESEARCH  OF  PROPERTIES  OF  SOLUTIONS  OF  CROSS-DIFFUSION 
MODEL  OF  REACTION  DIFFUSION  WITH  DOUBLE  NONLINEARITY 

N 

Consider  in  the  domain  Q={(t,x):  0<t<°o,  xe  R  [parabolic  system  of  cross-diflusion 


LLL  =  v(|x|"  |vMf  \p~2  V  u2‘ )  +  kx  (u,  -  <), 


du2 

~dt 


=v(l 


x\n  | VwJ  \P~2  V u"'2 )  +  jfcj  (w2  -  u22 ), 


(18) 


Ml|r=0  ’  ^2 1 1=0  ’ 

which  describes  process  of  reaction-diffusion  in  the  nonlinear  bipropellant  environment  which  diffusion 
coefficients  are  equal  |x|"  |Vw*  |P  2  >  |x|”  |Vm2  |P  2;  numeric  parameters  ml,m2,n,p,Pl,P2  -  positive  real  numbers, 

,  p2  >  0 ,  M,  =  Mj (t, x)  >  0  ,  u2=u2(t,x)>  0  -  density. 

In  this  work  it  is  investigated  properties  of  solutions  of  an  initial  value  problem  for  the  system  of  biological 
population  with  double  nonlinearity.  The  main  research  technique  is  self-similar  approach.  We  will  construct  a  self-similar 
set  of  equations  by  method  of  nonlinear  splitting  [17]. 

Substitution  in  (18) Ul(t,x)  —  e  U2(t,x)  —  e  kl'v2{r(t),  x)  will  lead  (18)  to  the  form: 


=  v(|x|"  |Vv* \P~2  Vv”1 )  - a^v* , 
=  V  ||x|"  |Vv*  |P_2  Vv'"2  j 


dr 
dv2 
.  dr 


(19) 


-a{Cblv22 


Vt  L=0=  Vio(jc)  ,  v2  ,=0=  v20(x). 


Here 


T(t)  = 


+( p—2)kki  —k{  ]r 


0[m2kl+(p-2)kk1-k2V 


mlk2  +  (p-2)kkl- k{  m2kl  +(p  —  2 )kk2  - k2 
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_  —  (p  — 2)kk1  -  m, k2  ^  _  k2/?,  —(p  —  2 )kkl  - nt ,k2 

1  n\k2  +  (p  —  2)kkt  —  k,  2  m2k{  +{p  —  2  )kkl  —  k 0 

Then  the  solution  of  system  (19)  is  sought  in  the  form 

Vl(t,x)  =  Vi(T)wl(T(t),<p(\x\)),  V2(t,x)  =V2{T)w2{T{t),(p(\x\j) 

yt(T)  =  (T0  +  T)~r'  ,  V2(T)  =  (T0  +  T)71  ,Tq>  0, 


1101 


b2+ 1 


Where  at  b,  =0,  b2  =0:  y,  = - ,  Y2  = - ,atb,  ^0,  b2  ^0:  Y,  =— 1 - ,  Y2  = 

12  '  A,-1  Pi- 1  1,2  a-i  A2-i 

Then  for  W;-  (T,  (/?(|a'|)),  i  =  1,2  get  the  system  of  equations: 


=  <P 


d 

r 

rQS-\ 

dwk 

p~2  dw 

dtp 

V 

V 

dtp 

dip 

d 

f 

ras~l 

dwk 

/,_2  m2 

ow  2 

dtp 

(p 

dip 

dip 

+  ^1(w,-w1A), 
+  Y2(w2-w^), 


(20) 


(21) 


where  y/{  =  - 


and  If/,  =  - 


Z)k  - 1]  -  Y2m !  )r,  (1  -  ^[(p  -  2)k-Y\  -  yng  )t2 

Here  at  p  >  n  :  (p(\x\)  =  \x\Pl  /  p{,  px=(p~n)l  p ,  s  =  pN  /  (p  -  n)  , 

and  at  p  =  n  :  <p( |x|)  =  ln(|x|)  , 

Self-similar  solution  of  system  (21)  looks  like 
W^Titltp)  =  /.(£),  £  =  tp(\x\)  I  ■ 

Then  substituting  (22)  in  (21)  relatively  f{  (p)  get  a  system  of  self-similar  equations 


dg 

dg 

where  If 


df k 


d% 

df2 


p-2 


d% 


p-2 


df2  ‘  )  +  +  n  (f  -f  ft)  =  0 

d^2df^h  h  }  ’ 

df"l\  %df2  sc  c d  \  a 
^2—) +  -^2. +  «,(/,  -  fft  )  =  0. 
d$  2  dig  2  2  2 

_ 1 _ 

{l-Yl[{p-2)k-\^r2ml)  ^  (l-y2[(p- 2)k-\]-ylm2) 

The  system  (23)  has  approximate  solution  of  a  look 

f\  =  A(a  -  £)+n' ,  J2  =  B(a  -  tg)^2 , 


1 


(22) 


(23) 
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where 

n  _(p-m(p-2)-nh\  n  _(p-m(p-2)-m2\ 

1  [k(p-2)f -mxm2  2  [k(p-2)f  -n\m2 

Theorem 

Let  Uj  (0,  x)  <  Zj  (0,  x\  X  G  ,  i  =  1,2  Then  for  the  solution  of  a  task  (18)  in  area  Q 

0 ,  X)  <  z{  (f ,  x)  =  (T  +  t)~n  /i  (£)  ,  U2  (t,  x)  <  Z2  (t,  X )  =  (T  +  0”^  /2  (^)  ,  £,  =  9(|x|)  /  T1/p  assessment  takes 
place. 

Here  the  functions  defined  above. 

In  case  n{  >0 >0,  g>0  applying  a  method  [17]  to  the  solution  of  the  equation  (23)  we  will  receive  the 
following  functions 

3(f)  =  (a-#)/\  ft(f)  =  (a-f  V, 

Theorem 

Limited  solution  of  a  system  (23)  at  ^  — >  a^has  asymptotics  /.(£)-#  . 

In  case  n{  >  0 ,n2  >0,  q  <  0  for  (23)  we  have 

*,(#)  =  (a +  £  )\  z2(f)=(a  +  f  r, 

where  <2  >  0  q  =  [k( p-  2)  —  l]2  —  mxm2 . 

Theorem 

At  ^  — >  +°°  task  solution  (23)  disappearing  on  infinity  the  solution  has  an  asymptotics  f.  (C,),  i  —  1, 2  . 

4.  CONCLUSIONS 

The  properties  stated  above  because  of  the  theorem  comparison  of  the  solution  established,  an  asymptotics  of  self¬ 
similar  solutions,  including  for  a  case  of  fast  diffusion  received.  Based  on  the  found  solutions,  numerical  calculations 
carried  out.  The  research  of  qualitative  properties  of  a  system  (18)  allowed,  to  execute  a  numerical  experiment  depending 
on  values,  the  logging-in  numerical  parameters.  For  this  purpose  as  an  initial  approximation  an  asymptotics  of  solutions 
used. 
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APPENDIX 


^  =  V  ( \x\n  |Vwf  r2  Vw2mi )  +  kx  (u,  -  uf' ) 
du ’  =  V  (|jc|"  | Vm*  I'"'  \u[”2  |  +  /Cj(m2  -  m2a  ) 


Fast  diffusion.  As  an  initial  approximation,  it  is  necessary  to  take: 


1 


1 


ux(x,t)  =  (T  +  T(t))  y'(a  +  ^ryi  ,u2(x,t)  =  (T  +  T(t))  Yl(a  +  %7)  " \yx  =——,y2  =— — - 

Pi  ~  *  Pl~^ 


7  =  — , £  =  <p(\x I)  /  Typ  p  >  n  :  (p(\x\)  =  | x\P'  / Pi,Pi=(p-n)/ p, ,h  = 
p-l  1  1 


<7 


(p-l)[k(p-2)-m2]  2 


,  q  =  [&(/?  —  2)]-  —  m,m2  ft,  >0,n2  >0,q<0 


Table  1 


Values  of  Parameters 


'.a  =f  *1.=  =0-972, 
=0-972 


4a  =2,  =1-046. 

=1-046 


^  =  3,^™,  =1,096, 
x2lim  =  1.096 


'^=4,^  =1.132, 
x,  =  1 .132 


<_  =  5!x1„=1.162, 
x,ot  =  1.!62 


Wj  =  5.9,  np  =7  p  =  3 
eps  =  10-3 
w,  =0.312>  0 
n,  =  0.365  >  0 
q  =  -41.05  <  0 

A  =  2,  A  =3 

A:  =  0.5n=0.3 


m,  =  8,  m2  =  7,  />  =  3 
e/w  =  10  “3 
ft,  =  0.255  >0 
n2  =  0.218  >0 
g  =  - 55  <  0 
A  =2  A  =3 
k-2 
n=0.5 


•  Fast  diffusion.  As  an  initial  approximation,  it  is  necessary  to  take 


7- 


ul(x,t)  =  (T  +  r(t ))  r,(a  +  i;r)n'  ,u2(x,t)  =  (T  +  r{t))  r2(a  +  ^y)  "2,yl  =— l—  ,72 

Pi  ~  *-  Pi  ~  *- 

-1  >  >  a  a 


q  =  [k{p  —  2)]2  -m{m2  n{  >0 ,n2  >0,q<0 
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Table  2 
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Values  of  Parameters 


;  =1.  *!„,.=  2-752, 
2.752 


^.  =  2,x,„„  =  3.152, 
=3.152 


=  3,  x,  =3.476, 


=  4,*1„„=3.742, 
=  3.742 


t  =5.  x  =3.974. 

tiax  Hinax 

=3.974 


=  10-3 
?2j  =  0255  >  0 

A  =0218  >0 

<?=  -55  <  0 
/?.  =2  fr=3k=2 
n=3 


TMj  =5.9,  m2=l  p=  3 
eps=  10~3 
A  =0.312  >0 
w3=  0.365  >0 
^  =  -41 .05  <  0 

A=2=  A=3 
A:  =0.5 

n=3 _ 


•  Slow  diffusion.  As  an  initial  approximation,  it  is  necessary  to  take: 


«,<.(,/)  =  (7  +  r(0)- " (a - ,  «,(x,r)  =  (T + r(0rK (fl ,  r,  = 


P-T  Yl  A-l 


7  =  — ,  £  =  p(|x|)  /  p>«:  p(|*|)  =  |-*f'  /  A ,  A  =  (p  -  n)  /  p , 


_  (p-l)[^(j7-2)-»7,]  _  (/?  — l)[fc(/7~  2)  — 777,] 


■  ,  72 


,  q  =  [k(/7  —  2)]2  —  tn fit 2  n\  >0,n2>0,q>0 . 


Table  3 


Values  of  Parameters 

m,  =0.5,  m2  =0.7,  y?=5 

=  10‘3 

«!  =  0.73  >  0 

n,  =  0.698  >0 
?  =  24.65  >0 
A  =2,  A  =3 

i  =  2i=0.5 

mi  =  2,  m:  =4,  £=10 

epj  =  10"3 

i\  =1242>0 
n,  =  0.65 
?  =30.44  >0 

A  =5,  A  =7 

k  =  0.9  n=0.5 

m,  =  02,  m.  =li,  £=7 
e/K  =  10"3 

n,  =0.654>0 
n2  =  0.558  >0 
^  =  80.7  >  0 
A  =2,  A  =3 
k  =  2 

n=2 


=3.  ^  =1.096. 
=1.096 


^=4=^  =  1.132, 
^  =  1.132 


•  Slow  diffusion.  As  an  initial  approximation,  it  is  necessary  to  take 
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ul(x,t)  =  (T  +  mrn(a-^)+n', 


u2(x,t)  =  (T  +  T(t)yn(a-er)*, 


1  I 

A7!'  72~A^' 


7  =  — D—  =  ^(Ixl)  /  r7,  p  =  n  :  (p{ |x|)  =  ln(|x|) 

p  —  1 


_  (j7-l)[fc(/7-2)-m,]  _  (^-l)[^(p-2)-m2] 


,72 


,  q  =  [&(/?  -  2)]2  -  mjm2  «j  >  0,n2  >  0,q  >  0. 


Table  4 


Values  of  Parameters 


ml  =0.5,  m2  =0.1,  p  =  5 
eps  =  10  "3 
«1  =  0.73>0 
n2  =  0.698  >0 
^  =  24.65  >0 
A=2,  A=3 

k  =  2  n=5 


mi  =2,  m2  =4,  p= 10 
=  10  "3 

n,  =1.242  >0 
n:  =0.65 
?  =  30 .44  >  0 
A  =5,  A  =7 
k  =  03  n-10 
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